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OH  Laser  Induced  Fluorescence  (LIF)  and  picosecond  (ps),  broadband  Coherent  Anti-Stokes  Raman  Spec¬ 
troscopy  (CARS)  are  used  for  time-resolved  temperature  and  time-resolved,  absolute  OH  number  density 
measurements  in  lean  H2-air,  CH4-air,  C2H4-air,  and  C3H8-air  mixtures  in  a  nanosecond  (ns)  pulse  dis¬ 
charge  cell/plasma  flow  reactor.  The  premixed  fuel-air  flow  in  the  reactor,  initially  at  T0  =  500  K  and 
P=  lOOtorr,  is  excited  by  a  repetitive  ns  pulse  discharge  in  a  plane-to-plane  geometry  (peak  voltage 
28  kV,  discharge  gap  10  mm,  estimated  pulse  energy  1.25  mj/pulse),  operated  in  burst  mode  at  10  kHz 
pulse  repetition  rate.  In  most  measurements,  burst  duration  is  limited  to  50  pulses,  to  preclude 
plasma-assisted  ignition.  The  discharge  uniformity  in  air  and  fuel-air  flows  is  verified  using  sub-ns-gated 
images  (employing  an  intensified  charge-coupled  device  camera).  Temperatures  measured  at  the  end  of 
the  discharge  burst  are  in  the  range  of  T  =  550-600  K,  using  both  OH  LIF  and  CARS,  and  remain  essentially 
unchanged  for  up  to  10  ms  after  the  burst.  Time-resolved  temperature  measured  by  CARS  during  plasma- 
assisted  ignition  of  H2-air  is  in  good  agreement  with  kinetic  model  predictions.  Based  on  CARS  measure¬ 
ment,  vibrational  nonequilibrium  is  not  a  significant  factor  at  the  present  conditions. 

Time-resolved,  absolute  OH  number  density,  measured  after  the  discharge  burst,  demonstrates  that  OH 
concentration  in  C2H4-air,  C3H8-air,  and  CH4  is  highest  in  lean  mixtures.  In  H2-air,  OH  concentration  is 
nearly  independent  of  the  equivalence  ratio.  In  C2H4-air  and  C3H8-air,  unlike  in  CH4-air  and  in  H2-air, 
transient  OH-concentration  overshoot  after  the  discharge  is  detected.  In  C2H4-air  and  C3H8-air,  OH  decays 
after  the  discharge  on  the  time  scale  of  ~0.02-0.1  ms,  suggesting  little  accumulation  during  the  burst  of 
pulses  repeated  at  10  kHz.  In  CH4-air  and  H2-air,  OH  concentration  decays  within  ~0.1-1.0  ms  and  0.5- 
1.0  ms,  respectively,  showing  that  it  may  accumulate  during  the  burst. 

The  experimental  results  are  compared  with  kinetic  modeling  calculations  using  plasma/fuel  chemistry 
model  employing  several  H2-air  and  hydrocarbon-air  chemistry  mechanisms.  Kinetic  mechanisms  for  H2- 
air,  CH4-air,  and  C2H4-air  developed  by  A.  Konnov  provide  the  best  overall  agreement  with  OH  measure¬ 
ments.  In  C3H8-air,  none  of  the  hydrocarbon  chemistry  mechanisms  agrees  well  with  the  data.  The  results 
show  the  need  for  development  of  an  accurate,  predictive  low-temperature  plasma  chemistry/fuel  chem¬ 
istry  kinetic  model  applicable  to  fuels  C3  and  higher. 

©  2013  The  Combustion  Institute.  Published  by  Elsevier  Inc.  All  rights  reserved. 


1.  Introduction 

Nonequilibrium  plasmas  formed  by  nanosecond  (ns)  duration 
pulsed  discharges  have  shown  significant  effect  on  ignition  delay 
reduction  and  flame  stabilization  [1-9].  The  peak  reduced  electric 
field,  defined  as  the  ratio  of  electric  field  and  total  gas  number  den¬ 
sity,  (E/N),  in  these  discharges  is  of  the  order  of  102-103Td 
(1  Td  =  10-17  V  cm2),  which  enhances  electronic  excitation  and 
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dissociation  of  gas  molecules,  resulting  in  efficient  generation  of 
radical  species.  The  radicals  readily  react  with  the  fuel  and  oxidizer 
molecules,  resulting  in  chain  branching  reactions  and  heat  genera¬ 
tion,  which  can  eventually  lead  to  ignition.  To  analyze  quantita¬ 
tively  the  impact  of  radicals  produced  in  the  plasma  on  fuel 
oxidation  and  ignition,  time-resolved  measurements  of  gas  tem¬ 
perature  and  radical  concentrations  are  required.  Recently,  laser 
diagnostics  has  been  applied  to  quantitative  analysis  of  plasma 
assisted  ignition/combustion,  such  as  OH  Planar  Laser  Induced 
Fluorescence  (PLIF)  [10,11],  saturated  OH  LIF  [12]  and  O  atom 
Two-photon  Absorption  LIF  (TALIF)  [13,14].  In  our  previous  work, 
various  laser  diagnostics  methods  have  been  used  for  temporally 
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and  spatially  resolved  measurements  in  ns  discharges  in  fuel-air 
mixtures,  including  absolute  NO,  OH,  and  O  atom  concentrations 
by  LIF  [15-18],  rotational  and  N2  vibrational  temperatures  by 
Coherent  Anti-Stokes  Raman  Spectroscopy  (CARS)  [19-21],  and 
rotational  temperature  by  OH  LIF  [18]. 

In  our  previous  measurements  of  ignition  delay  and  OH  concen¬ 
tration  in  H2-air  mixtures,  two  discharge  operation  regimes  have 
been  used.  In  the  first  approach,  a  ns  pulse  discharge  operating 
at  a  high  pulse  repetition  rate  was  maintained  through  the  ignition 
process  [17,22].  Kinetic  modeling  of  the  plasma  chemical  ignition 
process  demonstrated  that  in  this  case  ignition  time  had  a  fairly 
weak  dependence  on  the  details  of  kinetic  mechanism.  This  occurs 
due  to  continuous  generation  of  radicals  during  the  discharge, 
overlapping  with  ignition,  which  largely  “masks”  such  depen¬ 
dence.  In  the  second  approach,  the  discharge  burst  was  terminated 
shortly  before  the  onset  of  ignition  [18].  At  these  conditions,  elec¬ 
tronically  excited  molecules  and  atoms  generated  during  the  dis¬ 
charge  decay  rapidly,  such  that  oxidation  of  fuel-air  mixtures  is 
controlled  primarily  by  low-temperature  reactions  of  species  in 
their  ground  electronic  states,  with  the  pool  of  radicals  (primarily 
H,  O,  and  OH)  generated  during  the  discharge.  Therefore,  this  ap¬ 
proach  uncouples  plasma  chemical  reactions  from  “conventional” 
chemical  reactions,  making  analysis  of  kinetics  more  straightfor¬ 
ward.  In  this  case,  ignition  of  H2-air  mixtures,  initially  at 
T0  =  500  K,  occurred  when  the  temperature  at  the  end  of  the  dis¬ 
charge  burst  exceeded  a  threshold  of  approximately  T/^  700  K, 
about  200  K  lower  than  autoignition  temperature  predicted  by  ki¬ 
netic  modeling  [22].  Based  on  ignition  time,  time-resolved  temper¬ 
ature,  and  time-resolved  absolute  OH  concentration 
measurements  using  OH  LIF  [17,18,22],  insight  into  the  effect  of 
super-equilibrium  concentrations  of  radical  species  generated  in 
low-temperature  H2-air  plasmas  has  been  obtained. 

However,  quantitative  study  of  plasma  assisted  ignition  kinetics 
by  a  repetitive  ns  pulse  discharge  could  not  be  easily  extended  to 
hydrocarbon  fuels,  due  to  nonuniform  ignition  caused  even  by 
moderate  plasma  nonuniformity  [18,22].  In  the  present  work,  this 
problem  is  circumvented  by  reducing  the  number  of  pulses  in  the 


discharge  burst,  as  well  as  the  equivalence  ratio.  This  is  done  to 
study  the  mechanism  of  radical  generation  (with  OH  being  the  pri¬ 
mary  focus)  in  mildly  preheated  fuel-air  mixtures  excited  by  a 
repetitively  pulsed  ns  discharge,  as  well  OH  production  and  decay 
after  the  discharge,  before  ignition  would  occur.  The  objectives  of 
the  present  work  are:  (i)  to  use  both  OH  LIF  thermometry  and  pico¬ 
second  (ps)  CARS  for  time-resolved  rotational  temperature  mea¬ 
surements  of  fuel-air  mixtures  after  a  ns  pulse  discharge  burst, 
and  (ii)  to  measure  absolute,  time-resolved  OH  concentrations  in 
H2-air  and  hydrocarbon-air  mixtures  excited  by  the  discharge  after 
it  is  turned  off.  Thus,  plasma  chemical  reactions  and  “conven¬ 
tional”  chemical  reactions  are  separated  in  time,  similar  to  the  ap¬ 
proach  used  in  our  previous  work  [18,22].  These  measurements  are 
conducted  for  several  different  fuels,  over  a  range  of  equivalence 
ratios.  The  results  are  compared  with  kinetic  modeling  calcula¬ 
tions,  to  assess  the  applicability  of  H2-air  and  hydrocarbon-air 
chemistry  mechanisms  at  the  present  low-temperature  conditions. 

2.  Experimental 

The  experimental  setup  used  in  the  present  work,  shown  sche¬ 
matically  in  Fig.  1,  is  similar  to  the  one  used  in  our  previous  studies 
[18].  The  discharge  cell/plasma  flow  reactor  consists  of  a  280-mm- 
long,  22  mm  x  10  mm  rectangular  inner  cross  section  quartz  chan¬ 
nel  with  wall  thickness  of  1.75  mm.  Two  plane  quartz  windows, 
1  mm  thick  each,  are  fused  to  the  ends  of  the  channel,  providing 
optical  access  in  the  axial  direction.  A  60-mm-long,  right-angle 
fused  silica  prism  is  placed  along  the  channel  to  provide  optical  ac¬ 
cess  from  the  side  (see  Fig.  1).  The  entire  assembly  is  heated  in  a 
tube  furnace  (Thermcraft,  Ltd.,  with  152-mm-diameter,  254-mm- 
long  heated  section),  to  improve  plasma  stability.  A  1-m-long 
quartz  tube  coil  inlet  (see  Fig.  1)  preheats  the  fuel-air  flow  to  the 
furnace  temperature,  which  was  verified  by  thermocouple  mea¬ 
surements.  Two  6-mm-diameter  quartz-to-stainless-steel  adaptors 
connect  the  reactor  to  the  gas  delivery  system.  Fuel  and  air  flow 
rates  through  the  reactor  are  metered  by  MKS  mass  flow  control- 
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Fig.  1.  Schematic  diagram  of  the  preheated  plasma  assisted  ignition  test  cell  and  OH  LIF  apparatus. 
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lers  (<0.4%  error  at  the  experimental  flow  rates,  according  to  the 
manufacturer’s  specifications).  Fuel  and  air  flows  are  premixed  be¬ 
fore  entering  the  cell. 

Two  14  mm  x  60  mm  rectangular  plate  copper  electrodes, 
rounded  at  the  edges,  are  placed  on  the  top  and  bottom  of  the 
quartz  channel,  as  shown  in  Fig.  1,  and  held  in  place  by  ceramic 
clamps.  A  1.6-mm-thick  high-temperature  perfluoroelastomer 
sheet  (Kalrez,  DuPont)  is  placed  between  each  electrode  and  the 
channel  wall,  to  reduce  air  gaps  and  prevent  corona  discharge  out¬ 
side  the  cell.  Thus,  the  total  distance  between  the  electrodes  is 
approximately  17  mm,  with  the  discharge  gap  (distance  between 
top  and  bottom  quartz  walls)  of  10  mm.  The  electrodes  are  con¬ 
nected  to  an  FID  GmbH  FPG  60-100MC4  pulse  generator  (peak 
voltage  up  to  30  kV,  pulse  duration  of  5  ns,  repetition  rate  up  to 
100  kHz).  In  the  present  work,  the  pulser  is  operated  in  repetitive 
burst  mode,  producing  bursts  of  50  pulses  at  a  pulse  repetition  rate 
of  10  kHz  and  burst  repetition  rate  of  5  Hz. 

A  schematic  diagram  of  the  OH  LIF  setup  is  also  shown  in  Fig.  1. 
Briefly,  the  532  nm  output  from  an  Nd:YAG  laser  (Continuum, 
Model  Powerlite  8010)  is  used  to  pump  a  tunable  dye  laser  (Laser 
Analytical  Systems,  Model  LDL  20505),  the  output  of  which  is  fre¬ 
quency  doubled  by  a  p  Barium  Borate  (BBO)  crystal  (wavelengths 
shown  in  Fig.  1  correspond  to  OH  A-X  (1,0)  pumping).  While  the 
flash  lamp  of  the  Nd:YAG  laser  is  operated  at  10  Hz,  Q-switching 
controlling  the  laser  output  is  maintained  at  5  Hz,  to  match  the 
burst  repetition  rate  of  the  pulser  and  to  probe  OH  after  each  dis¬ 
charge  burst.  A  delay  generator  is  used  to  control  the  delay  be¬ 
tween  the  discharge  burst  and  the  laser  pulse.  A  half  wave  plate 
and  thin  film  polarizer  pair  is  used  to  control  the  UV  beam  energy 
for  LIF  operation  in  the  linear  regime  (below  5  pj/pulse).  The  UV 
beam  energy  is  monitored  using  a  photodiode  placed  after  a  beam 
splitter  (see  Fig.  1). 

OH  concentration  and  rotational  temperature  (assumed  to  be  in 
equilibrium  with  the  translational  temperature)  are  measured  on 
the  channel  centerline  by  OH  LIF  in  lean  CH4-,  C2H4-,  C3H8-,  and 
H2-air  mixtures,  each  at  four  different  equivalence  ratios,  at  flow 
velocity  through  the  reactor  of  u  «  40  cm/s.  The  purity  of  the  fuels 
is  as  follows:  99.0%  for  CH4  (<0.02%  of  02,  H20,  and  C02  combined, 
<0.2%  of  N2),  99.5%  for  C2H4  and  C3H8,  and  99.999%  for  H2.  The 
compressed  air  cylinder  used  is  Ultra  Zero  grade,  with  19.5- 
23.5%  of  02  and  <2  ppm  H20.  A  low  discharge  burst  repetition  rate, 
5  Hz,  is  chosen  to  keep  the  delay  between  the  bursts  longer  than 
the  residence  time  of  the  flow  in  the  plasma,  ~0.1  s,  such  that 


the  flow  in  the  reactor  would  experience  only  one  burst.  The  low 
flow  velocity  also  improves  flow  preheating  in  the  inlet  coil  and  re¬ 
duces  pressure  drop  across  the  reactor. 

For  relative  OH  concentration  measurements,  the  Q.i(3)  line  in 
the  OH  A-X  (1,0)  band  is  used  to  maximize  signal-to-noise  ratio 
at  temperatures  below  600  K.  For  OH  thermometry,  four  additional 
lines  in  the  (1,0)  band  system  are  used  to  generate  a  more  accurate 
Boltzmann  fit  for  rotational  temperature,  including  the  Pi(l ),  Qi(2), 
Q.i(4),  and  Q.i(5)  lines.  The  choice  of  the  excitation  lines,  as  well  as 
LIF  signal  processing  was  discussed  in  detail  in  Ref.  [17].  Additional 
details  of  LIF  diagnostics  are  discussed  in  Section  4.  LIF  is  collected 
from  the  center  of  the  discharge,  a  region  approximately  15-mm- 
long,  via  a  right  angle  prism  and  fused  silica  lenses,  and  directed 
to  a  photomultiplier  tube  (PMT)  (Hamamatsu,  Model  R106)  in 
combination  with  a  small  monochromator  (Horiba,  Model  H-10), 
which  is  used  to  collect  fluorescence  signal  only  from  (0,0)  and 
(1,1)  transitions.  The  LIF  signal  is  integrated  over  a  100-ns  gate 
and  analyzed  in  real  time  during  the  experiment  using  a  program¬ 
mable  digital  oscilloscope  with  1-GHz  bandwidth  (LeCroy,  Model 
WAVEPRO  7100A). 

Relative  OH  LIF  concentrations  measured  using  (1,0)  excitation 
are  calibrated  using  the  LIF  signal  from  Ri(4)  excitation  in  the  (0,0) 
band  and  Rayleigh  scattering  at  the  wavelength  close  to  the  center 
of  the  (0,0)  fluorescence  band,  308  nm,  as  discussed  in  detail  in 
Section  4.  For  calibration  using  Rayleigh  scattering,  the  flow  chan¬ 
nel/discharge  cell  has  been  modified  to  accommodate  the  entrance 
and  exit  windows  set  at  Brewster’s  angle.  Note  that  using  Rayleigh 
scattering  calibration  requires  no  change  in  optical  alignment, 
therefore  mitigating  uncertainties  encountered  with  the  use  of  an 
atmospheric-pressure  Hencken  burner  flame  [18]. 

Rotational  temperature  and  vibrational  temperature  of  N2  on 
the  centerline  of  the  discharge  cell/plasma  flow  reactor  are  also 
measured  by  ps  Coherent  Anti-Stokes  Raman  Spectroscopy  (CARS) 
diagnostics  described  in  detail  in  Ref.  [21]  and  shown  schemati¬ 
cally  in  Fig.  2.  Briefly,  in  this  system  a  broadband  Stokes  beam  is 
generated  by  an  in-house  fabricated  modeless  dye  laser,  with  con¬ 
version  efficiency  of  approximately  6%.  The  dye  laser  is  pumped  by 
an  Ekspla  SL-333  Nd:YAG  laser,  generating  output  pulses  ~150  ps 
in  duration,  with  energy  output  of  up  to  120  mj/pulse  at  532  nm. 
The  Nd:YAG  laser  also  generates  the  pump/probe  beams  for  the 
CARS  mixing.  The  dye  laser  output  is  centered  near  604  nm,  with 
a  full  width  at  half  maximum  of  approximately  5-6  nm.  In  the 
present  work,  Unstable-resonator  Spatially  Enhanced  Detection 


Pump  /  Probe 
beam 


Fig.  2.  Schematic  diagram  of  ps  USED-CARS  apparatus. 
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Table  1 

Dominant  primary  radical  species  generation  processes  in  the  plasma  [12,15,18,21]. 


Process 

Rate  (cm3/s) 

Al 

N2  +  e“  =  N2(A3Z,  B3n,  C3n,  a'1!])  +  e“ 

<ra 

A2 

N2  +  e“  =  N(4S)  +  N(4S)  +  e“ 

G 

A3 

02  +  e“  =  0(3P)  +  O^P^D)  +  e“ 

G 

A4 

N2(c3n)  +  o2  =  n2  (B3n)  +  o2 

3.0  X  10-10 

A5 

N2(a/1Z)  +  02  =  N2  (B3n)  +  02 

2.8  x  lO”11 

A6 

N2(B3n)  +  02  =  N2  (A3£)  +  02 

3.0  x  10“10 

A7 

N2(A3£)  +  02  =  N2  +  O  +  O 

2.5  x  10“12 

HI 

H2  +  e-  =  H  +  H  +  e- 

<72 

H2 

N2(a/1Z)  +  H2  =  N2  +  H  +  H 

2.6  x  10“n 

H2 

N2(B3n)  +  H2  =  N2(A3£)  +  h2 

2.5  x  10“n 

H4 

N2(A3£)  +  H2  =  N2  +  H  +  H 

4.4  x  10-1°  exp(-3170/T) 

H5 

0(1D)  +  H2  =  H  +  OH 

1.1  x  lO”10 

Ml 

CH4  +  e“  =  CH3  +  H  +  e“ 

G 

M2 

N2(A3£)  +  CH4  =  N2  +  CH3  +  H 

1.2  X  10-1°  exp(-3500/T) 

M3 

N2(B3n)  +  ch4  =  n2  +  ch3  +  h 

3.0  x  lO”10 

M4 

N2(c3n)  +  ch4  =  n2  +  ch3  +  h 

5.0  x  10“10 

M5 

N^a'1^)  +  CH4  =  N2  +  CH3  +  H 

3.0  x  10“10 

El 

C2H4  +  e“  =  products3  +  e“ 

G 

E2 

N2(A3£)  +  C2H4  =  N2  +  C2H3  +  H 

9.7  x  10“n 

E3 

N2(B3n)  +  C2H4  =  n2  +  C2H3  +  H 

3.0  x  10“10 

E4 

N2(C3n)  +  c2h4  =  n2  +  C2H3  +  H 

3.0  x  lO”10 

E5 

N^a'1^)  +  C2H4  =  N2  +  C2H3  +  H 

4.0  x  10“10 

PI 

N2(A3£)  +  C3H8  =  N2  +  C3H6  +  H2 

1.2  x  10“12 

P2 

N2(B3I1)  +  c3h8  =  n2  +  c3h6  +  h2 

3.0  x  10“10 

P3 

N2(c3n)  +  c3h8  =  n2  +  c3h6  +  h2 

3.0  x  lO”10 

P4 

N2(a/1£)  +  C3H8  =  N2  +  C3H6  +  H2 

4.0  x  10“10 

a  Calculated  by  the  Boltzmann  solver  from  the  experimental  cross  sections. 
b  Sum  of  electronic  excitation  cross  sections  (b3Z,  tTz,  c3n,  a3Z,  Un,  and  d3n). 
c  Three  dissociation  channels  into  C2H3  +  H,  C2H2  +  H2,  and  C2H2  +  H  +  H. 


(USED)  CARS  phase  matching  has  been  employed,  with  longitudi¬ 
nal  spatial  resolution  of  ~5  mm.  For  this  geometry,  the  532  nm 
pump/probe  beam  is  enlarged  and  the  center  portion  is  blocked, 
creating  an  annulus,  as  shown  in  Fig.  2.  This  beam  is  then  com¬ 
bined  coaxially  with  the  Stokes  beam  by  a  dichroic  mirror  and  fo¬ 
cused  into  the  test  section.  After  re-collimation,  dichroic  mirrors 
isolate  the  CARS  signal  beam,  which  is  then  directed  into  a  0.75- 
m  spectrometer  (see  Fig.  2).  At  the  spectrometer  exit  plane,  a 
2.3x  relay  lens  magnification  system  is  used  to  image  the  dispersed 
CARS  signal  onto  an  electron-multiplying  charge  coupled  device 
(CCD)  camera  (Andor  Newton).  The  combination  of  the  2.3x  magni¬ 
fication  and  3600  line/mm  grating  results  in  a  spectral  resolution 
of  ~0.4  cm-1,  which  is  sufficient  to  partially  resolve  the  rotational 
structure  in  the  Q-branch  of  N2.  Details  of  data  reduction,  as  well  as 
rotational  and  vibrational  temperature  inference  from  the  CARS 
spectra  are  discussed  in  Ref.  [21]. 

3.  Plasma  chemistry  model  and  nanosecond  pulse  discharge 
model 

To  obtain  insight  into  kinetics  of  plasma/chemical  fuel  oxida¬ 
tion  and  ignition,  we  use  a  kinetic  model  developed  in  our  previous 
work  [15-20,22].  Briefly,  the  model  incorporates  nonequilibrium 
air  plasma  chemistry  [23],  expanded  to  include  hydrocarbons 
and  hydrogen  dissociation  processes  in  the  plasma  [16,24].  In  the 
present  work,  the  plasma  chemical  reaction  list  was  also  expanded 
to  incorporate  C3H8  dissociation  reactions  by  electronically  excited 
N2  [25].  The  dominant  radical  species  generation  processes  in  H2- 
air  and  hydrocarbon-air  and  plasmas  are  listed  in  Table  1.  The  plas¬ 
ma  chemistry  model  is  coupled  with  one  of  two  “conventional” 
hydrogen-oxygen  chemistry  mechanisms  developed  by  Popov 
[24]  and  Konnov  [25],  and  one  of  the  three  hydrocarbon  chemistry 
mechanism,  GRI  Mech  3.0  [23],  USC/Wang  mechanism  [24],  or 
Konnov  mechanism  [29].  Note  that  the  mechanisms  have  been 
developed  and  validated  for  relatively  high-temperatures, 
significantly  higher  than  found  in  the  present  experiments,  and 


may  well  be  inaccurate  at  the  present  conditions.  Assessing  appli¬ 
cability  of  these  mechanisms,  used  as  a  starting  point  for  develop¬ 
ment  of  a  nonequilibrium  plasma  assisted  combustion  mechanism, 
is  one  of  the  objectives  of  the  present  work. 

The  list  of  air  plasma  chemistry  processes  and  hydrogen-oxy¬ 
gen  chemical  reactions  incorporated  in  the  present  model  is  given 
in  Ref.  [16]  and  Refs.  [24,26],  respectively.  The  species  concentra¬ 
tion  equations  are  coupled  with  the  two-term  expansion  Boltz¬ 
mann  equation  for  plasma  electrons.  The  model  incorporates  the 
energy  equation  for  the  temperature  on  the  discharge  centerline 
[22]  as  well  as  quasi-one-dimensional  flow  equations,  with  heat 
transfer  to  the  walls  (dominant  energy  loss  mechanism  at  the  pres¬ 
ent  conditions)  included.  The  model  is  validated  using  measure¬ 
ments  of  O-atom  concentration  (TALIF)  [15,16],  temperature 
(rotational  CARS)  [19,20],  OH  concentration  (LIF)  and  ignition  time 
[17,18,22]  in  ns  discharges  in  air,  and  mixtures  of  air  and  CH4, 
C2H4,  and  H2. 

Key  parameters  controlling  plasma  chemistry  in  the  ns  pulse  dis¬ 
charge  include  the  reduced  electric  field,  F/N,  and  coupled  pulse  en¬ 
ergy.  Energy  coupled  to  the  plasma  using  two  different  ns  pulse 
generators,  including  the  FID  pulser  used  in  the  present  work,  was 
measured  in  air  in  a  low-temperature  cell  outside  of  the  furnace, 
over  a  wide  range  of  pressures  [30].  It  has  been  shown  that  at 
P  =  1 0-1 00  torr,  coupled  pulse  energy  is  nearly  independent  of  pulse 
repetition  rate  (at  v  =  1  -40  kHz),  remains  nearly  constant  during  the 
pulse  burst  (up  to  ~1 00  pulses),  and  is  proportional  to  discharge 
pressure,  i.e.,  energy  coupled  per  molecule  remains  constant.  These 
experimental  results  are  in  good  agreement  with  the  analytic  model 
of  energy  coupling  in  a  ns  pulse  discharge  [31  ],  which  incorporates 
key  effects  of  pulsed  breakdown  and  sheath  development  on  the 
ns  time  scale.  The  model  predicts  pulse  energy  coupled  to  the  plasma 
vs.  pulse  voltage  waveform,  discharge  geometry,  pressure,  and  tem¬ 
perature.  Time-dependent  reduced  electric  field  in  the  plasma  dur¬ 
ing  the  discharge  pulse,  predicted  by  the  model  of  Ref.  [31],  is 
fairly  complex,  with  F/N  rapidly  dropping  during  breakdown  [30], 
which  introduces  an  uncertainty  in  the  “effective”  F/N  value  at 
which  energy  is  coupled  to  the  plasma.  In  the  present  plasma  chem¬ 
istry  model,  F/N  in  the  plasma  is  approximated  as  a  Gaussian  pulse 
with  the  same  time  constant  as  the  experimental  voltage  waveform, 
and  the  peak  corresponding  to  the  breakdown  value  predicted  by  the 
ns  discharge  model  [31  ],  approximately  (F/N)peak  =  600  Td,  with  the 
coupled  pulse  energy  of  1.25  mj/pulse  (0.46  meV/molecule).  At  the 
present  conditions,  the  model  predictions  are  not  very  sensitive  to 
the  effective  F/N  value.  Reducing  peak  F/N  by  a  factor  of  3,  to  200 
Td,  while  keeping  the  pulse  energy  the  same,  results  in  a  fairly  mod¬ 
est  increase  of  predicted  OH  number  density  after  the  burst,  by 
approximately  20%.  The  uncertainty  of  the  predicted  coupled  energy 
value  is  approximately  ±50%,  due  to  a  significant  uncertainty  in  the 
dielectric  constant  of  high-temperature  dielectric  layers  (Kalrez).  In 
the  experiment,  the  uncertainty  is  primarily  due  to  the  effect  of  cor¬ 
ona  discharges  forming  near  the  electrodes  outside  the  test  cell, 
which  reduces  the  energy  coupled  to  the  plasma  inside  the  cell. 
The  value  of  1.25  mj/pulse,  used  in  the  present  work,  provides  the 
best  agreement  with  the  data.  Varying  the  pulse  energy  within  the 
estimated  uncertainty  limit  results  in  predicted  [OH]  variation  by 
approximately  ±20%. 

It  is  necessary  to  point  out  that  the  plasma  chemistry  model 
used  in  this  work  also  incorporates  a  set  of  NOx  formation  and  de¬ 
cay  reactions  (list  of  these  reactions  for  air  plasma  is  given  in  Ref. 
[16]).  NO  may  affect  radical  species  chemistry  at  low  temperatures, 
primarily  via  recombination  reaction  NO  +  O  (+M)  ->  N02  (+M), 
which  reduces  O  atom  number  density.  NO  may  also  affect  OH 
number  density,  e.g.,  via  a  rapid  reaction  H02  +  NO  ^  N02  +  OH. 
At  the  present  conditions,  the  predicted  NO  mole  fraction  after 
the  burst  is  of  the  order  of  ~150ppm,  significantly  higher 
compared  to  OH  mole  fraction,  ~10ppm.  However,  varying  NO 
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fraction  in  the  initial  fuel-air  mixture  in  the  range  of  0-300  ppm 
resulted  in  a  weak  change  of  the  predicted  OH  number  density, 
by  1-2%.  Therefore  at  the  present  low-temperature  conditions, 
NO  chemistry  is  not  expected  to  play  an  important  role  in  fuel  oxi¬ 
dation  kinetics. 


4.  OH  LIF  and  LIF  thermometry 

4.1.  Basics 


Operating  in  linear  excitation  regime,  the  time-integrated  fluo¬ 
rescence  signal,  5/,  can  be  expressed  as  follows: 


Sf(v)  =  n0HfBbuE(/)jg(v)l—fi. 


(1) 


where  n0H  is  OH  number  density  (cm-3);  fB  is  the  Boltzmann  factor 
for  OH  molecules  in  the  absorbing  state,  calculated  using  the  rota¬ 
tional  energy  term  values  for  X2n  v=0  state  from  Ref.  [32]; 
b\2  =  B12IC  is  the  Einstein  absorption  coefficient  (cm2  J-1  cm-1),  ta¬ 
ken  from  LIFBASE  [33];  E  is  the  laser  energy  (J);  </>;  is  the  fluores¬ 
cence  quantum  efficiency  (or  quantum  yield);  g(v)  is  the 
normalized  overlap  integral  (1/cm-1),  i.e.,  a  convolution  of  laser 
lineshape  and  absorption  lineshape;  /  is  the  path  length  of  the  laser 
sampled  by  the  collection  optics  (cm);  Q  is  the  solid  angle  of  detec¬ 
tion  (sr);  and  p  represents  the  overall  conversion  coefficient,  which 
consists  of  the  efficiency  of  the  collection  optics  and  response  of  the 
detection  system.  For  (1,0)  excitation,  the  quantum  yield  can  be  de¬ 
rived  by  integrating  the  rate  equation  over  the  entire  LIF  pulse, 
while  also  taking  into  account  spontaneous  emission  (Ai  and  A0), 
electronic  quenching  (Ch  and  Qo),  and  vibrational  energy  transfer 
(VET)  from  v7  =  1  to  v7  =  0  (V10): 


E\A\  +  £qAqUio/ (Aq  +  Qq) 

V10  +A\  +  Q.! 


(2) 


In  Eq.  (2),  £0  and  £1  are  signal  collection  efficiencies  for  emission 
from  v7  =  0  and  v7  =  1  of  the  A2Z+  state,  respectively,  which  account 
for  differences  in  fluorescence  trapping  and  spectral  sensitivity  of 
the  detection  system  for  different  transitions.  It  should  be  noted 
that  VET  acts  to  redistribute  the  A-state  population  (from  v7  =  1 
to  v7  =  0),  and  the  equation  can  be  simplified  with  the  following 
two  approximations  that  eliminate  the  dependence  on  the  VET 
rate:  (1 )  A0  +  Qo  «  A\  +  Q.i  (which  is  not  unreasonable  based  on  lit¬ 
erature  quenching  rates,  see  Refs.  [34,35])  and  (2)  £oA0  «  s^A  1  (also 
not  unreasonable,  see  Ref.  [33]).  Thus, 


 £0^0 


(3) 


In  this  work,  the  quantum  yield  is  calculated  from  the  measured 
fluorescence  decay  time,  t  =  1/(A0  +  Qo).  We  note  that  the  quantum 
yield  estimate  for  (1,0)  excitation  can  be  compared  with  measure¬ 
ments  using  A-X(0,0)  excitation. 

Without  knowing  the  laser  linewidth  accurately,  calculation  of 
the  normalized  overlap  integral  in  Eq.  (1 ),  g(  v\  may  result  in  signif¬ 
icant  uncertainty.  In  the  present  work,  performing  the  excitation 
scan  over  an  isolated  transition,  and  integrating  the  fluorescence 
signal  over  a  wide  spectral  range,  g{  v)  becomes  unity,  and  therefore 
no  knowledge  of  laser  linewidth  is  required  to  reduce  the  data.  In 
this  case,  Eq.  (1)  becomes 

Sf  =  noufsbuE^jl^P-  (4) 

For  OH  LIF  thermometry,  (1,0)  excitation  is  preferred  for  better 
rejection  of  laser  scattering.  Rotational  temperature  is  inferred 
from  the  slope  of  the  Boltzmann  fit  of  fluorescence  intensities  from 
the  five  excitation  transitions  listed  in  Section  2. 


4.2.  Calibration 

In  order  to  obtain  a  calibration  factor  for  the  LIF  data  taken 
using  (1,0)  excitation,  experiments  employing  two  of  the  four 
equivalence  ratios  for  each  fuel-air  mixture  listed  in  Section  2 
were  repeated  using  (0,0)  excitation.  These  experiments  have  been 
conducted  in  the  modified  test  cell  with  Brewster  windows  (see 
Section  2).  Since  VET  is  not  involved  in  (0,0)  excitation  (assuming 
no  upward  VET  from  v7  =  0  to  v7  =  1 ),  the  fluorescence  quantum 
yield,  calculated  from  the  measured  fluorescence  decay  time  using 
Eq.  (3),  is  more  accurate.  The  rest  of  parameters  in  Eq.  (4),  neces¬ 
sary  to  calculate  the  absolute  OH  concentration,  i.e.,  the  product 
of  l,  Q,  and  /?,  are  determined  from  Rayleigh  scattering,  as  follows: 

W-b©”' '<NE>'  (5> 

In  Eq.  (5),  h  is  the  Planck  constant  (J  s);  v  is  the  laser  wavelength 
(cm-1);  c  is  the  speed  of  light  (cm/s);  (|g)  is  the  Rayleigh  scattering 
differential  cross  section  (cm2/sr),  calculated  using  dispersion  rela¬ 
tions  for  depolarization  from  Ref.  [36]  and  for  index  of  refraction 
from  Ref.  [37];  N  =  P/kT  is  the  gas  number  density  (cm-3);  and  sR 
is  the  collection  efficiency  of  the  Rayleigh  scattering.  By  varying 
the  pressure  of  a  buffer  gas  (N2  or  air)  inside  the  same  test  cell  as 
used  for  OH  LIF  experiments,  and/or  the  laser  pulse  energy,  a  linear 
plot  of  SRayieigh  vs.  (NE)  was  obtained,  the  slope  of  which  was  used  to 
calculate  the  collection  constant,  IQ/3,  in  Eq.  (5).  Then,  absolute  OH 
concentrations  were  calculated  for  the  (0,0)  excitation  cases.  Final¬ 
ly,  by  comparing  signal  levels  to  those  with  known  OH  concentra¬ 
tion,  the  rest  of  relative  OH  concentrations  measured  using  (1,0) 
excitation  were  put  on  an  absolute  scale.  This  calibration  scheme 
is  discussed  in  greater  detail  in  Ref.  [38].  Based  on  a  standard  error 
propagation  analysis,  and  by  taking  into  account  uncertainties  in 
measured  laser  energy,  inferred  temperature,  shot-to-shot  signal 
variation,  and  uncertainties  in  the  cross  sections  of  various  pro¬ 
cesses  involved  (absorption,  emission,  quenching,  vibrational  en¬ 
ergy  transfer,  and  Rayleigh  scattering),  the  combined  uncertainty 
in  the  OH  number  density  is  estimated  to  be  ±20%. 


5.  Results  and  discussion 

5.1.  Discharge  uniformity 

Plasma  uniformity  is  difficult  to  quantify  but  critical  for  isolat¬ 
ing  nonequilibrium  chemistry  effects  from  thermal  heating  occur¬ 
ring  in  filaments  and  hot  spots.  The  approach  used  for  sustaining 
diffuse,  nearly  zero-dimensional  ns  pulse  discharges  in  a  plane- 
to-plane  geometry  was  discussed  in  detail  in  our  previous  work 
[18,22].  Before  quantitative  OH  LIF  measurements  were  performed, 
plasma  uniformity  was  verified  qualitatively  by  taking  images 
employing  an  intensified  charge-coupled  device  (ICCD)  camera 
during  individual  discharge  pulses.  Time-resolved,  broadband,  sin¬ 
gle-shot  images  taken  during  a  high  voltage  ns  pulse  discharge  in 
air  at  T0  =  500  K  and  P  =  50  torr,  using  a  PI-MAX  ICCD  camera  and 
a  UV  lens  (Sodern  100-mm  f/2.8),  are  shown  in  Fig.  3.  The  camera 
gate  is  set  at  490  ps  and  synchronized  with  the  high  voltage  pulse 
generator  to  open  during  the  last  pulse  in  a  10-pulse  discharge 
burst,  with  t  =  0  representing  the  beginning  of  the  pulse.  The  tim¬ 
ing  of  the  gate  is  indicated  on  the  top  of  each  image.  The  images 
capture  primarily  nitrogen  second  positive  band  emission,  N2 
C3n  ->  B3n.  From  the  images  in  Fig.  3,  it  is  evident  that  the  plasma 
is  generated  in  the  entire  volume  between  the  electrodes  and  re¬ 
mains  diffuse  and  nearly  uniform  until  it  decays  completely.  Total 
emission  intensity  (a  sum  of  intensity  counts  over  all  pixels  in  each 
ICCD  image)  decays  very  rapidly  after  each  discharge  pulse,  on  the 
time  scale  of  a  few  ns,  demonstrating  that  excited  radiating  species 
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Fig.  3.  Broadband  single-shot  ICCD  images  of  the  plasma  during  a  single  discharge  pulse  (pulse  #10  in  the  burst)  in  air.  To  =  500K,  p  =  50  torr,  v  =  10kHz.  Camera 
gate  =  490  ps,  timing  of  the  gate  is  shown  on  top  of  each  image.  Intensity  of  images  taken  at  t  <  10  ns  is  enhanced  for  illustrative  purposes. 


(mainly  N2  C3I1  state)  are  quenched  rapidly.  Figure  4  shows  ICCD 
images  taken  at  the  same  conditions  as  in  Fig.  3,  through  the  win¬ 
dow  at  the  end  of  the  flow  channel.  To  isolate  plasma  emission 
from  the  reflections  from  the  channel  walls,  images  were  taken 
with  and  without  the  discharge,  and  a  rectangular  region  corre¬ 
sponding  to  the  channel  cross  section  (22  mm  x  10  mm)  was  ex¬ 
tracted  from  the  original  images,  as  shown  in  Fig.  4.  Again,  it  can 
be  seen  that  the  plasma  is  diffuse  and  occupies  the  entire  volume 
between  the  electrodes. 

As  discussed  in  our  previous  work  [22],  increasing  pressure  and/ 
or  adding  H2  or  hydrocarbons  to  the  air  plasma  makes  the  dis¬ 
charge  less  uniform  and  results  in  formation  of  well-defined  fila¬ 
ments,  due  to  development  of  ionization  -  heating  instability. 
This  effect  is  especially  well  pronounced  in  hydrocarbon-air  mix¬ 
tures  at  pressures  of  ~100  torr  and  above,  although  mild  preheat¬ 
ing  to  T0  =  400-500  K  reduces  the  nonuniformity  considerably. 
Figure  5  shows  ICCD  images  (50-ns  exposure)  of  discharges  in 
air,  H2-air  (</>  =  0.3),  and  C2H4-air  (</>  =  0.3)  preheated  to 
T0  =  500  K,  taken  through  the  window  in  the  end  of  the  channel. 
These  images  illustrate  that  the  plasma  remains  diffuse  and  fila¬ 
ment-free  at  these  conditions.  As  in  Fig.  4,  the  images  show  the  en¬ 
tire  channel  cross  section  (22  mm  x  10  mm).  Note,  however,  that 


ns  pulse  plasma-assisted  ignition  in  CH4-air  and  C2H4-air  mixtures 
still  exhibits  a  well-defined  propagating  flame  front  [22],  such  that 
near-homogeneous  ignition  could  not  be  achieved.  In  the  present 
work,  this  problem  is  circumvented  by  employing  short  discharge 
bursts  (50  pulses  at  the  pulse  repetition  rate  of  v  =  10  kHz,  far  less 
than  needed  for  ignition)  and  operating  in  fuel-lean  mixtures,  with 
the  main  emphasis  of  measuring  OH  number  density  in  the  plasma 
prior  to  ignition.  Plasma  uniformity  at  these  conditions  has  been 
verified  for  three  representative  fuel-air  mixtures  (CH4-,  C2H4-, 
and  H2-air  mixtures  at  </>  =  0.3),  at  three  different  pressures  (100, 
300  and  500  torr),  as  shown  in  Fig.  6.  As  discussed  in  Section  2, 
LIF  signal  is  collected  from  the  region  in  the  center  of  the  plasma 
approximately  1.5  cm  long,  about  %  of  the  length  of  the  discharge. 
This  eliminates  the  effect  of  discharge  nonuniformity  that  may  be 
found  near  electrode  edges  at  high  pressures  [22]. 

5.2.  OH  LIF  thermometry  and  comparison  to  CARS  temperature 
measurements 

A  typical  excitation  spectrum  of  the  five  OH  A-X  (1,0)  transi¬ 
tions  used  in  the  present  work  is  shown  in  Fig.  7.  The  spectrum 
is  taken  2  ps  after  the  end  of  the  last  pulse  in  an  H2-air  mixture 
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Fig.  4.  Top:  single  shot  ICCD  images  of  the  flow  channel  without  the  discharge  and  with  the  discharge,  taken  through  the  window  at  the  end  of  the  channel.  Extracted  region 
corresponds  to  the  channel  cross  section  (22  mm  x  10  mm).  Bottom:  images  of  the  plasma  during  a  single  discharge  pulse  (pulse  #10  in  the  burst)  in  air  at  the  conditions  of 
Fig.  3,  taken  through  the  same  window  and  showing  the  entire  channel  cross  section.  Camera  gate  490  ps,  timing  of  the  gate  is  shown  on  top  of  each  image. 


1600 


Z.  Yin  et  al.  /  Combustion  and  Flame  7  60  (2013)  1594-1608 


Fig.  5.  Single  shot  ICCD  images  of  the  plasma  in  and  fuel-air  mixtures  preheated  to  T0  =  500  K,  taken  through  the  window  in  the  end  of  the  channel,  v  =  10  kHz,  camera 
gate  =  50  ns,  images  show  the  entire  channel  cross  section  (22  mm  x  10  mm). 


Fig.  6.  Single  shot  ICCD  images  of  the  plasma  in  fuel-air  mixtures  illustrating 
discharge  uniformity.  T0  =  500  K,  0  =  0.3,  v  =  10  kHz,  pulse  #10.  Camera  gate  =  50  ns. 


at  (j)  =  0.12,  T0  =  500  K,  and  P  =  100  torr,  excited  by  a  v  =  10  kHz, 
50  pulse  discharge  burst.  This  is  close  to  minimum  time  delay  for 
which  LIF  signal  is  not  affected  by  EMI  from  the  discharge.  In 
Fig.  7,  each  experimental  data  point  represents  a  median  of  60  laser 
shots  taken  at  each  scanning  step.  As  discussed  in  our  previous 
work  [18],  Voigt  line  shape  profiles  are  used  to  fit  the  experimental 
spectra  to  isolate  overlapping  spectral  lines,  which  appear  in  all  Q- 
branch  transitions.  After  isolating  the  main-branch  transitions  and 
integrating  them  over  the  same  wavelength  range,  symmetric  to 
the  centerline  wavelength,  rotational  temperature  is  inferred  from 
the  Boltzmann  plot,  such  as  shown  in  Fig.  8.  Each  of  the  excitation 
scans  shown  in  Fig.  7  is  repeated  three  times  to  improve  accuracy. 
Only  the  major  transitions  are  included  in  the  Boltzmann  fit.  At 
some  conditions,  including  the  one  shown  in  Fig.  8,  data  from 
the  Pi(l)  and  Q.i(2)  transitions  deviate  somewhat  from  the  straight 


line  formed  by  Qi(3),  Q.i(4)  and  Q.i(5)  transitions.  If  Pi(l)  and  Qi(2) 
are  removed  from  Fig.  8,  the  best  fit  rotational  temperature  would 
be  Trot  =  557  ±  25  K,  which  differs  by  30  K  from  the  best  fit  temper¬ 
ature  obtained  using  a  5  rotational  line  fit  (Trot  =  587  ±  27  K).  This 
relatively  small  discrepancy  could  be  due  to  P^l)  line  being  rela¬ 
tively  weak,  as  well  as  a  strong  overlap  between  Q.i(2)  with 
Q2i(2)  lines,  which  may  introduce  uncertainty  in  data  processing 
and  Voigt  line  fitting.  Of  course,  some  fitting  error  may  be  incurred 
with  run  to  run  variations  (since  there  may  in  fact  be  some  temper¬ 
ature  variation  from  run  to  run).  Nonetheless,  this  difference  in 
best  fit  rotational  temperature  values  is  close  to  the  statistical 
uncertainty  (95%  confidence  interval)  of  the  Boltzmann  fit  used 
for  temperature  inference.  In  most  cases,  transitions  involving 
low  rotational  levels  (e.g.,  Pi(l)  and  Qi(2))  are  needed  for  better 
rotational  temperature  sensitivity. 

Figure  9a  shows  a  100-shot  accumulation  N2  ps  CARS  spectrum 
taken  in  air  at  P=  100  torr  and  To  =  500I<  (without  the  plasma), 
along  with  the  best-fit  synthetic  spectrum  generated  using  the 
Sandia  National  Laboratories  CARSFIT  program  [39].  The  confi¬ 
dence  interval  of  temperature  inference  was  determined  by  taking 
100  CARS  spectra  at  these  well-characterized  conditions,  100-shot 
accumulation  each.  For  each  of  these  spectra,  rotational  tempera¬ 
ture  was  inferred  using  a  least  squares  fit  synthetic  spectrum  cal¬ 
culated  by  CARSFIT.  The  best  fit  rotational  temperature,  T  =  486  K, 
standard  deviation,  and  95%  confidence  interval,  ~15K  at  these 
reference  conditions,  were  obtained  from  statistical  analysis  of 
these  100  temperature  values.  The  relative  uncertainty  in 
rotational  temperature  is  expected  to  improve  with  increasing 
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Fig.  7.  Typical  OH  LIF  excitation  spectrum  used  for  temperature  measurements.  Spectrum  taken  2  ps  after  the  last  pulse  in  the  burst  (v  =  10  kHz,  50  pulses).  H2-air,  T0  =  500  K, 
P=100  torr,  0  =  0.12. 
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Fig.  8.  Boltzmann  plot  obtained  from  the  OH  LIF  excitation  spectrum  shown  in 
Fig.  7  (H2-air,  T0  =  500I<,  P=  100  torr,  0  =  0.12,  v  =  10kHz,  50  pulses).  Best  fit 
rotational  temperature  is  T  =  587  ±  27  K. 
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Fig.  9.  (a)  100-Shot  accumulation  ps  CARS  spectrum  in  air  at  P=100torr  and 
T0  =  500  K,  Tfit  =  486  K,  95%  confidence  interval  ~15  K  and  (b)  100-shot  accumula¬ 
tion  ps  CARS  spectrum  in  0  =  0.4  H2  -  air  plasma  at  P  =  92  torr  during  ignition, 
Tfit  =  1250  I<. 


temperature  due  to  increasing  populations  of  higher  rotational  lev¬ 
els,  which  are  more  readily  resolved  in  the  Q-branch  spectra.  Fig¬ 
ure  9b  plots  a  100-shot  accumulation  N2  ps  CARS  spectrum  taken 
in  4>  =  0.4  H2-air  plasma  at  P  =  92  torr  during  ignition,  and  the  best 
fit  synthetic  spectrum  with  Tnt  =  1250  K.  In  Fig.  9b,  N2(p=  1)  vibra¬ 
tional  band  is  also  detected  at  Raman  shifts  below  2300  cm-1.  The 
two  experimental  spectra  shown  in  Fig.  9  represent  two  extremes 
of  the  temperature  observed  in  the  present  experiments,  500  K  and 
1250 1<  (during  ignition).  As  can  be  seen,  in  both  cases,  the 


synthetic  spectra  are  in  good  agreement  with  the  experimental 
spectra.  100-shot  accumulation  CARS  spectra  were  used  to  obtain 
each  temperature  value  discussed  below. 

Figure  10  compares  time-resolved  rotational  temperature  mea¬ 
sured  by  OH  LIF  in  a  C2H4-air  mixture  at  4>  =  0.09,  T0  =  500  K, 
P  =  100  torr,  after  a  v  =  10  kHz,  50-pulse  discharge  burst  with  rota¬ 
tional  temperature  measured  by  ps  CARS  at  the  same  conditions. 
Temperatures  measured  by  these  two  techniques  are  close  to  each 
other,  showing  that  the  temperature  in  the  discharge  increases  by 
AT  =  50-80  K  (from  the  initial  value  of  T0  =  500  K)  and  remains 
nearly  constant  for  up  to  1  ms  after  the  discharge  burst.  In  fact, 
additional  temperature  measurements  for  very  long  delays  after 
a  50-pulse  burst,  plotted  in  Fig.  12a  and  discussed  below,  demon¬ 
strated  that  the  temperature  remains  approximately  the  same  for 
up  to  ~10  ms  after  the  burst.  These  results  are  also  consistent  with 
the  temperature  rise  in  the  discharge  predicted  by  the  kinetic  mod¬ 
el  at  these  conditions  (equaling  AT  ~  100  K),  also  shown  in  Fig.  10. 
A  relatively  modest  temperature  rise  in  the  discharge,  as  well  as  its 
weak  dependence  on  the  equivalence  ratio  are  also  observed  in 
other  fuel-air  mixtures  at  the  same  initial  temperature,  pressure, 
and  discharge  burst  duration,  as  shown  in  Fig.  1 1.  The  temperature 
rise  during  the  burst,  averaged  over  different  fuels  and  equivalence 
ratios,  is  ATavg  «  70  K.  Again,  the  temperatures  derived  from  OH 
LIF  data  (triangles,  Fig.  11)  are  consistent  with  the  CARS  results 
(squares)  and  fuel-air  plasma  chemistry  model  predictions  (lines). 
Based  on  these  results,  the  same  temperature  value,  averaged  over 
different  fuels  and  equivalence  ratios,  Tavg  =  570  K,  was  used  to  cal¬ 
culate  the  Boltzmann  factor  to  infer  absolute,  time-resolved  OH 
number  densities  from  the  LIF  data. 

Figure  12a  plots  temperatures  measured  by  CARS  after  the  dis¬ 
charge  burst  in  H2-air  at  4>  =  0.4,  T0  =  500  K,  P  =  92  torr,  and 
v  =  10  kHz,  for  three  different  burst  durations  (50-120  pulses).  Also 
plotted  in  Fig.  12a  are  the  model  predictions  at  these  conditions 
(solid  lines),  which  are  in  very  good  agreement  with  CARS  data 
up  to  fairly  long  delay  times  after  the  burst,  20  ms,  although  the 
model  again  systematically  overpredicts  the  temperature  at  the 
end  of  the  discharge  burst,  by  about  50  K.  From  the  results  shown 
in  Fig.  12a,  it  is  apparent  that  threshold  temperature  at  the  end  of 
the  discharge  burst,  above  which  ignition  occurs,  is  700  K, 
approximately  200  K  below  the  normal  autoignition  temperature 
[22].  This  result  confirms  our  previous  temperature  measurements 
by  N2  emission  spectroscopy  [22],  although  the  results  of  Ref.  [22] 
were  not  spatially  resolved  and  had  more  uncertainty.  Below  this 
threshold,  ignition  does  not  occur  and  the  temperature  gradually 
decays  to  pre-excitation  baseline  value  of  T0  =  500  K.  In  addition, 
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Fig.  10.  Comparison  of  temperatures  measured  by  OH  LIF  and  by  CARS  after  a 
v  =  10  kHz,  50-pulse  discharge  burst  in  a  C2H4-air  mixture  at  T0  =  500  K,  P  =  100  torr, 
0  =  0.09. 
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Fig.  11.  Comparison  of  temperatures  measured  by  OH  LIF  and  by  CARS  in  a 
v  =  10  kHz,  50-pulse  discharge  burst  (2  ps  after  the  last  pulse)  in  different  fuel-air 
mixtures  at  T0  =  500  K,  P=  100  torr.  Model  predictions  are  shown  as  solid  lines. 
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Fig.  12.  Comparison  of  ps  CARS  temperature  measurements  after  the  discharge 
burst  in  a  ^  =  0.4  H2-air  mixture  at  ro  =  500K,  P  =  92  torr  and  v  =  10kHz  with 
kinetic  modeling  calculations,  for  different  burst  durations  (50-120  pulses).  Plot  (b) 
is  a  semi-log  graph  (for  the  time  axis)  of  plot  (a)  with  the  addition  of  OH  number 
density  predicted  by  the  model.  Model  predictions  are  shown  as  lines.  Solid  lines, 
temperature;  dashed  lines,  OH  number  density. 


Fig.  12b  plots  the  same  temperatures  vs.  time  as  a  semi-log  graph, 
along  with  OH  concentrations  predicted  by  the  model  (shown  by 
dashed  lines).  It  can  be  seen  that  below  ignition  threshold,  OH  gen¬ 
erated  during  the  discharge  burst  also  gradually  decays  and  does 
not  exhibit  an  overshoot  typical  for  ignition. 

Although  CARS  measurements  demonstrated  that  fairly  signifi¬ 
cant  vibrational  excitation  of  nitrogen  occurs  in  a  v  =  10  kHz  ns 
pulse  burst  discharge  sustained  in  the  present  cell  in  room  temper¬ 
ature  air  [22],  N2  vibrational  temperature  measured  in  preheated 
air  (T0  =  500  K)  was  quite  low,  TV(N2)  «  850  K  for  50  pulses  and  Tv(- 
N2) «  1050  K  for  100  pulses,  as  shown  in  Fig.  13.  In  fuel-air  mix¬ 
tures,  N2  vibrational  temperature  after  a  50-pulse  burst  was  even 
lower,  near  the  detection  limit,  TV(N2) «  600-700  K.  This  demon¬ 
strates  that  at  the  present  experimental  conditions,  vibrational 
nonequilibrium  of  excited  fuel-air  mixtures  is  not  a  significant  fac¬ 
tor.  Since  estimated  characteristic  decay  time  of  excited  electronic 
states  of  nitrogen  at  the  present  conditions  is  quite  short,  t  ~  1  ps 
for  N2(A3Z),  and  even  shorter  for  N2(B3n),  N2(C3n)  and  N2(a7lZ) 
states,  chemical  kinetics  of  fuel-air  mixtures  on  longer  time  scales 
is  dominated  by  “conventional”  chemical  reactions  among  radical 
species  generated  in  the  discharge.  This  lends  the  present  experi¬ 
ment  to  studies  of  radical  reaction  kinetics  at  relatively  low  tem¬ 
peratures,  ranging  from  room  temperature  to  ignition  temperature. 

5.3.  OH  number  density  after  the  discharge 

Figure  14  plots  absolute,  time-resolved  [OH]  measured  in  lean 
H2-air  mixtures  at  T0  =  500  K  and  P  =  100  torr  vs.  time  delay  after 
a  v  =  10  kHz,  50-pulse  discharge  burst,  at  different  equivalence  ra¬ 
tios  ((f)  =  0.015-0.12).  Rotational  temperatures  measured  by  LIF 
thermometry  at  the  end  of  the  burst  (time  delay  2  ps  after  the 
end  of  the  last  pulse  in  the  discharge  burst)  at  these  conditions 
are  also  shown  in  Fig.  14  (circled).  It  can  be  seen  that  peak  [OH] 
after  the  burst,  (3.0-3.6)1013  cm-3,  is  not  sensitive  to  the  equiva¬ 
lence  ratio  in  this  range.  Also,  there  is  essentially  no  transient  [OH] 
rise  after  the  burst.  The  fairly  long  [OH]  decay  time,  ~0.5-1.0  ms, 
shows  that  OH  accumulates  during  the  10  kHz  burst  (tine  interval 
between  the  pulses  of  0.1  ms),  consistent  with  our  previous  results 
in  H2-air  at  discharge  pulse  repetition  rate  of  v  =  401<Hz  [13].  As 
discussed  above,  temperature  at  the  end  of  the  discharge  burst  is 
not  sensitive  to  the  equivalence  ratio  and  delay  time  after  the  dis¬ 
charge  (see  Figs.  10  and  11)  and  remains  within  the  range  of 
T=  540-580  K  (see  Fig.  14). 

Figure  14  also  compares  kinetic  model  predictions  with  the 
experimental  results,  using  two  different  “conventional”  H2-air 
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Fig.  13.  Vibrational  temperature  of  nitrogen  measured  by  CARS  after  a  v  =  10  kHz 
discharge  burst  in  air  T0  =  500  K,  P=  100  torr  (50  and  100  pulses). 
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Fig.  14.  Comparison  of  experimental  (y-axis  on  the  left)  and  predicted  (y-axis  on  the  right)  time-resolved,  absolute  OH  number  densities  after  a  v  =  10  kHz,  50-pulse 
discharge  burst  in  H2-air  mixtures  at  T0  =  500  K,  P  =  100  torr,  and  different  equivalence  ratios.  Data  points  for  temperatures  measured  at  the  end  of  the  burst  (time  delay  2  ps) 
are  circled,  (a)  Popov  mechanism,  (b)  Konnov  mechanism  and  (c)  Konnov  mechanism  with  modified  rate  of  reaction  OH  +  O  ^  H  +  02. 


chemistry  mechanisms  developed  by  Popov  [24]  (Fig.  14a)  and 
Konnov  [26]  (Fig.  14b).  Note  that  the  model  predictions  for  [OH] 
are  plotted  using  the  vertical  axis  on  the  right-hand  side  of  the  fig¬ 
ures,  while  the  experimental  data  are  shown  using  the  vertical  axis 
on  the  left-hand  side.  As  can  be  seen,  the  Popov  mechanism 


matches  weak  dependence  of  [OH]  on  the  equivalence  ratio  (by 
50-60%)  but  consistently  overpredicts  absolute  [OH]  and  some¬ 
what  overpredicts  the  [OH]  decay  time  at  the  higher  equivalence 
ratio  of  0  =  0.12  (see  Fig.  14a).  On  the  other  hand,  the  Konnov 
mechanism  matches  absolute  [OH]  concentration  as  well  as  decay 
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Fig.  15.  Comparison  of  experimental  (y-axis  on  the  left)  and  predicted  (y-axis  on  the  right)  time-resolved,  absolute  OH  number  densities  after  a  v  =  10  kHz,  50-pulse 
discharge  burst  in  CH4-air  mixtures  at  T0  =  500  K,  P=  100  torr,  and  different  equivalence  ratios.  Data  points  for  temperatures  measured  at  the  end  of  the  burst  (time  delay 
2  ps)  are  circled,  (a)  GRI  3.0  mechanism,  (b)  USC/Wang  mechanism  and  (c)  Konnov  mechanism. 


time  (except  at  0  =  0.12)  but  shows  a  somewhat  greater  depen¬ 
dence  of  peak  [OH]  on  the  equivalence  ratio  than  seen  in  the  exper¬ 
imental  results.  Kinetic  sensitivity  analysis  has  shown  that  the 
difference  between  the  two  mechanisms  is  almost  entirely  due  to 
the  difference  in  the  rate  coefficients  of  reaction  O  +  OH  ->  H  +  02, 


kr=  1.2  x  10-11  cm3/s  [19]  and  kr  =  5.3  x  10”11  cm3/s  [26]  at 
T  =  500  K  (note  that  the  former  value  violates  the  detailed  balance, 
since  the  rate  coefficients  used  for  the  reaction  H  +  02  -►  O  +  OH  in 
Refs.  [20,26]  are  close,  kf=  5.3  x  10“17cm3/s  and  /</=  4.9  x  10“17 
cm3/s,  respectively).  The  use  of  the  rate  coefficient  recommended 
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Fig.  16.  Comparison  of  experimental  and  predicted  time-resolved,  absolute  OH  number  densities  after  a  v  =  10kHz,  50-pulse  discharge  burst  in  C2H4-air  mixtures  at 
r0  =  500  K,  P=  100  torr,  and  different  equivalence  ratios.  Data  points  for  temperatures  measured  at  the  end  of  the  burst  (time  delay  2  ps)  are  circled,  (a)  GRI  3.0  mechanism, 
(b)  USC/Wang  mechanism  and  (c)  Konnov  mechanism. 


in  the  recent  review  [40],  kr  =  3.5  x  10-11  cm3/s  {kf=  3.3  x  10“17 
cm3/s)  in  the  Konnov  mechanism  results  in  [OH]  overprediction 
by  approximately  25%,  with  other  changes  being  insignificant 
(see  Fig.  14c). 


The  results  in  lean  CH4-air  mixtures  are  shown  in  Fig.  15.  Unlike 
in  H2-air  mixtures,  peak  [OH]  after  the  burst  is  reduced  as  the 
equivalence  ratio  increased,  by  about  50%  between  </>  =  0.03  and 
</>  =  0.24.  Similar  to  H2-air  mixtures,  almost  no  transient  [OH] 
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Fig.  17.  Comparison  of  experimental  and  predicted  time-resolved,  absolute  OH  number  densities  after  a  v  =  10kHz,  50-pulse  discharge  burst  in  C3H8-air  mixtures  at 
T0  =  500  K,  P=  100  torr,  and  different  equivalence  ratios.  Data  points  for  temperatures  measured  at  the  end  of  the  burst  (time  delay  2  ps)  are  circled,  (a)  GRI  3.0  mechanism, 
(b)  USC/Wang  mechanism  and  (c)  Konnov  mechanism. 


increase  is  observed,  except  for  a  modest  rise  at  </>  =  0.24.  The  [OH] 
decay  time  in  CH4-air  is  somewhat  shorter  than  in  H2-air,  ranging 
from  ~1  ms  at  </>  =  0.03  to  ~0.2  ms  at  </>  =  0.24.  Again,  this  suggests 
fairly  significant  OH  accumulation  during  the  discharge  burst  at 
the  pulse  repetition  rate  of  v  =  10  kHz.  At  the  higher  equivalence 
ratio,  4>  =  0.24,  a  modest  transient  [OH]  overshoot  after  the  burst 


is  apparent.  [OH]  after  the  discharge  burst  predicted  by  two  “con¬ 
ventional”  hydrocarbon-air  chemistry  mechanisms,  GRI  3.0  [27] 
and  USC/Wang  [28],  are  relatively  close  to  each  other  (see 
Fig.  15a  and  b),  both  exceeding  peak  measured  [OH]  by  60-100%. 
Both  mechanisms  reproduce  weak  transient  [OH]  rise  in  higher 
equivalence  ratio  mixtures,  with  Wang  mechanism  showing  better 
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agreement  with  [OH]  decay  after  the  burst.  The  Konnov  mecha¬ 
nism,  on  the  other  hand,  matches  absolute  [OH],  the  weak  tran¬ 
sient  rise,  and  the  [OH]  decay  relatively  well  (see  Fig.  15c). 
Comparison  with  the  modeling  calculations  shows  clearly  that 
both  peak  [OH]  and  decay  time  reduction  at  higher  equivalence  ra¬ 
tios  are  primarily  due  to  increase  in  net  rate  of  radical  removal  via 
fuel  oxidation  reactions,  since  coupled  discharge  pulse  energy  and 
the  net  rate  of  radical  generation  are  nearly  independent  of  the 
equivalence  ratio. 

An  additional  series  of  calculations  was  done  to  evaluate  the  po¬ 
tential  effect  of  hydrocarbon  impurities  in  99.0%  pure  methane. 
The  presence  of  trace  hydrocarbon  species,  such  as  such  as  ethyl¬ 
ene,  acetylene,  and  ethane,  rapidly  reacting  with  the  radicals  gen¬ 
erated  in  the  plasma  (primarily  with  O  and  H  atoms)  may  affect 
kinetics  of  plasma-assisted  methane  oxidation,  since  low-temper¬ 
ature  rates  for  CH4  reactions  with  0  and  H  atoms  are  fairly  slow. 
The  upper  bound  effect  of  these  species  on  low-temperature  meth¬ 
ane  oxidation  kinetics  was  estimated  using  kinetic  modeling  calcu¬ 
lations,  by  adding  1%  of  C2H2,  C2H4,  C2H6,  or  CO  to  CH4.  The  results 
demonstrate  that  C2H2,  C2H4,  or  C2H6  addition  results  in  less  than 
0.5%  difference  in  predicted  peak  temperature,  less  than  1%  differ¬ 
ence  in  predicted  quasi-steady-state  OH  number  density  at  the  end 
of  the  burst,  and  approximately  5%  difference  (reduction)  in  pre¬ 
dicted  OH  decay  time  after  the  burst.  As  expected,  OH  decay  time 
reduction  is  primarily  due  to  significantly  more  rapid  reaction 
rates  of  radicals  with  these  hydrocarbon  additives.  Adding  1%  CO 
did  not  result  in  a  detectable  change  of  the  predicted  temperature 
or  OH  number  density  at  the  present  conditions.  Although  higher 
methane  purity  is  indeed  desired  in  the  future  experiments,  we 
conclude  that  the  effect  of  hydrocarbon  impurities  (at  or  below 
1%  in  concentration)  in  methane  on  low-temperature  methane  oxi¬ 
dation  kinetics  is  fairly  minor. 

Figure  16  plots  the  results  in  lean  C2H4-air  mixtures.  It  can  be 
seen  that  in  the  entire  range  of  equivalence  ratios,  [OH]  after  the 
burst  first  increases,  by  up  to  60%,  before  eventually  decaying.  Sim¬ 
ilar  to  the  results  in  CH4-air  mixtures,  increasing  the  equivalence  ra¬ 
tio  reduces  peak  transient  [OH]  after  the  burst  by  approximately  a 
factor  of  2,  from  x  1013  cm-3  at  </>  =  0.05  to  x  1013  cm-3  at 
</>  =  0.36  (see  Fig.  16).  Finally,  [OH]  decay  time  is  much  shorter  than 
in  H2-air  and  CH4-air  and  is  reduced  considerably  as  the  equivalence 
ratio  is  increased,  from  ~0.1  ms  at  </>  =  0.05  to  ~0.03  ms  at  </>  =  0.36. 
This  shows  that  no  significant  OH  accumulation  occurs  during  the 
discharge  burst,  since  [OH]  decay  time  is  comparable  with  time  de¬ 
lay  between  the  pulses  in  a  10  kHz  burst.  Again,  these  trends  (both 
peak  [OH]  and  decay  time  reduction)  are  due  to  more  rapid  net  rate 
of  radical  removal  via  fuel  oxidation  reactions  at  the  higher  equiva¬ 
lence  ratio.  As  can  be  seen,  GRI  3.0  mechanism  overpredicts  both 
peak  [OH]  and  decay  time,  by  a  factor  of  3-4  (see  Fig.  16a).  The 
USC/Wang  mechanism  performs  noticeably  better,  although  peak 
[OH]  and  decay  time  are  still  underpredicted  by  about  a  factor  of  2 
(except  at  higher  equivalence  ratio,  </>  =  0.36,  where  the  model 
underpredicts  the  decay  time  by  about  a  factor  of  4,  see  Fig.  1 6b).  Fi¬ 
nally,  the  Konnov  mechanism  reproduces  both  peak  [OH]  and  its  de¬ 
cay  in  very  lean  mixtures  well,  although  the  agreement  is  not  as  good 
as  the  equivalence  ratio  is  increased  (see  Fig.  16c).  In  this  case,  peak 
[OH]  at  higher  equivalence  ratios  (</>  =  0.18  and  </>  =  0.36)  is  under¬ 
predicted  because  the  net  rate  of  C2H4  oxidation  process,  initiated 
by  C2H4  +  O  reactions,  is  overpredicted. 

Based  on  the  comparison  of  the  present  experimental  results 
with  kinetic  modeling,  the  Konnov  mechanisms  for  H2-air,  CH4- 
air,  and  C2H4-air  show  better  overall  agreement  with  the  data, 
compared  to  the  Popov,  GRI  3.0,  and  USC/Wang  mechanisms.  Note 
that  kinetic  modeling  shows  that  significant  transient  [OH]  rise 
after  the  burst,  observed  in  C2H4-air,  is  not  due  to  chain  branching 
processes,  which  remain  insignificant  at  these  low  temperatures. 
In  fact,  the  transient  rise  in  [OH]  is  produced  primarily  by  O  and 


H  atoms  generated  by  electron  impact  and  by  quenching  of  elec¬ 
tronically  excited  nitrogen  during  and  immediately  after  the  last 
pulse  in  the  burst  (within  ~1  ps  [14]). 

The  experimental  results  in  C3H8-air  are  very  similar  to  those  in 
C2H4-air  (compare  Figs.  16  and  17).  In  both  fuel-air  mixtures,  peak 
[OH]  is  reduced  as  the  equivalence  ratio  is  increased,  exhibiting  a 
transient  rise  after  the  burst,  and  the  [OH]  decay  rates  in  these 
two  fuel-air  mixtures  are  similar.  However,  kinetic  modeling  sig¬ 
nificantly  overpredicts  peak  [OH]  (by  a  factor  of  3  for  GRI  3.0,  a  fac¬ 
tor  of  6  for  USC/Wang,  and  a  factor  of  2  for  Konnov).  Also,  [OH] 
decay  time  is  overpredicted  considerably  by  all  three  mechanisms, 
and  the  significant  transient  rise  in  the  lean  mixture  (</>  =  0.04)  is 
not  reproduced  by  any  of  them.  In  particular,  predictions  from 
the  Konnov  mechanism,  which  performs  relatively  well  in  H2-air, 
CH4-air,  and  C2H4-air,  are  at  variance  with  the  experimental  results 
in  C3H8-air. 

We  believe  that  these  differences  are  primarily  due  to  uncer¬ 
tainties  in  “conventional”  C3  chemical  reaction  pathways  and  rates. 
Indeed,  the  fuel  fraction  in  the  leanest  C3H8-air  mixture  tested  (</> 
=  0.04)  is  only  0.18%  (compared  to  21%  of  air),  suggesting  that  dom¬ 
inant  plasma  chemical  radical  generation  processes  occur  during 
electron  impact  excitation  and  quenching  of  excited  N2  and  02, 
within  ~1  ps  after  the  burst.  Kinetics  of  these  processes  appear 
to  be  captured  by  the  present  plasma  chemical  model  (using  Kon¬ 
nov  mechanism)  relatively  well,  since  the  model  predictions  in 
very  lean  H2-air,  CH4-air,  and  C2H4-air  are  close  to  the  experiment 
values.  Thus,  the  present  data  show  the  need  for  development  of  an 
accurate,  predictive  plasma  chemistry/fuel  chemistry  kinetic  mod¬ 
el  that  would  be  applicable  to  fuels  C3  and  higher.  Also,  tracing  the 
cause  of  the  differences  between  different  fuel-air  kinetic  mecha¬ 
nisms  to  individual  species  and  reactions  requires  detailed  sensi¬ 
tivity  analysis,  which  will  be  undertaken  in  our  future  work. 


6.  Summary 

In  the  present  work,  calibrated  OH  Laser  Induced  Fluorescence 
(LIF)  is  used  for  time-resolved  temperature  and  time-resolved, 
absolute  OH  number  density  measurements  in  fuel-air  mixtures 
in  a  nanosecond  (ns)  pulse  discharge  cell/plasma  flow  reactor.  Rel¬ 
ative  OH  concentrations  are  put  on  an  absolute  scale  by  calibrating 
the  optical  collection  constant  using  Rayleigh  scattering.  Rotational 
temperature  is  inferred  from  OH  LIF  excitation  spectra  using  5  dif¬ 
ferent  excitation  transitions.  Rotational  temperature  and  N2  vibra¬ 
tional  temperature  in  the  reactor  have  also  been  measured  using 
broadband,  picosecond  Coherent  Anti-Stokes  Raman  Spectroscopy 
(CARS)  diagnostics.  The  premixed  fuel-air  flow  in  the  reactor  is  ex¬ 
cited  by  a  repetitive  ns  pulse  discharge  in  a  plane-to-plane  geom¬ 
etry,  operated  at  a  high  pulse  repetition  rate,  v  =  101<Hz.  The 
experiments  are  conducted  in  lean,  slowly  flowing  H2-air,  CH4- 
air,  C2H4-air,  and  C3H8-air  mixtures  preheated  to  T0  =  500  K  in  a 
tube  furnace,  at  a  pressure  of  P=  lOOtorr.  In  most  of  the  present 
measurements,  the  discharge  burst  duration  was  limited  to  50 
pulses,  insufficient  to  produce  plasma-assisted  ignition.  The  dis¬ 
charge  uniformity  in  air  and  fuel-air  flows  has  been  verified  using 
sub-ns-gated  images  from  an  intensified  camera  demonstrating 
that  a  diffuse,  volume  filling,  uniform  plasma  is  sustained  in  the 
entire  range  of  operating  conditions,  at  pressures  of  up  to 
P  =  500  torr. 

The  results  demonstrate  that  rotational  temperatures  inferred 
from  OH  LIF  spectra  and  partially  rotationally  resolved  CARS  spectra 
are  in  good  agreement,  providing  additional  validation  to  OH  LIF 
thermometry  diagnostics  used  in  the  present  paper,  as  well  as  in 
our  previous  work  [13].  Temperatures  at  the  end  of  the  discharge 
burst,  measured  by  these  two  techniques,  are  in  the  range  of 
T  =  550-580  K,  indicating  a  modest  temperature  rise  of  50-80  K 
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occurs  in  the  discharge  burst,  consistent  also  with  kinetic  modeling 
predictions.  The  temperature  after  the  burst  remains  essentially 
unchanged,  for  up  to  10  ms  delay.  N2  vibrational  temperatures  in¬ 
ferred  from  CARS  measurements  are  quite  low,  Tv  «  850  K  in  air 
and  Tv  «  600-700  K  in  fuel-air  after  a  50-pulse  burst,  demonstrat¬ 
ing  that  vibrational  nonequilibrium  is  not  a  significant  factor  at 
the  present  conditions. 

The  measured,  time-resolved,  absolute  OH  number  density 
demonstrates  that  OH  generated  in  C2H4-air  and  C3H8-air  mixtures 
is  highest  in  very  lean  mixtures  and  is  reduced  by  about  a  factor  of 
2  as  the  equivalence  ratio  is  increased  by  a  factor  of  ~8,  from  </> 
=  0.05  to  0.36  for  C2H4-air  mixture  and  from  </>  =  0.04  to  0.32  for 
C3H8-air  mixture.  In  these  two  fuels,  transient  OH  overshoot  after 
the  discharge  burst,  by  up  to  a  factor  of  2,  is  detected.  No  signifi¬ 
cant  overshoot  is  observed  in  CH4-air  and  in  H2-air  mixtures,  how¬ 
ever.  In  CH4-air  mixtures,  OH  concentration  is  also  reduced  as  the 
equivalence  ratio  is  reduced,  but  in  H2-air  mixtures,  OH  concentra¬ 
tion  after  the  burst  is  nearly  independent  of  the  equivalence  ratio. 
OH  decay  rate  after  the  burst  in  C2H4-air  and  C3H8-air  occurs  on  the 
time  scale  of  ~0.02-0.1  ms,  suggesting  little  accumulation  during 
the  burst  of  pulses.  In  CH4-air  and  H2-air,  however,  the  OH  decays 
within  ~0.1-1.0  ms  and  ~0.5-1.0  ms,  respectively,  suggesting  sig¬ 
nificant  OH  accumulation  in  lean  mixtures  during  the  burst. 

The  experimental  results  are  compared  with  kinetic  modeling 
calculations  using  a  plasma/fuel  chemistry  model  employing  sev¬ 
eral  different  “conventional”  H2-air  and  hydrocarbon-air  chemistry 
mechanisms.  Based  on  this  comparison,  kinetic  mechanisms  for 
H2-air,  CH4-air,  and  C2H4-air  developed  by  A.  Konnov  show  better 
overall  agreement  with  the  experimental  results,  compared  to 
H2-air  mechanism  developed  by  Popov,  as  well  as  the  GRI  3.0 
and  USC/Wang  mechanisms.  In  C3H8-air,  however,  none  of  the 
hydrocarbon  chemistry  mechanisms  agrees  well  with  the  data. 
The  present  results  show  the  need  for  development  of  an  accurate, 
predictive  low-temperature  plasma  chemistry/fuel  chemistry  ki¬ 
netic  model  applicable  to  fuels  C3  and  higher. 
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